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Abstract 

The concept of the nonholonomic deformation formulated recently for the AKNS family is ex- 
tended to the Kaup-Newell class. Applying this construction we discover a novel mixed integrable 
£f) ' hierarchy related to the deformed derivative nonlinear Schrodinger (DNLS) equation and found the 

exact soliton solutions exhibiting unusual accelerating motion for both its field and the perturbing 
functions. Extending the idea of deformation the integrable perturbation of the gauge related Chen- 
Lee-Liu DNLS equation is constructed together with its soliton solution. We show that, the recently 
proposed Lenells-Fokas (LF) equation falls in the deformed DNLS hierarchy, sharing the accelerating 
soliton and other unusual features. Higher order integrable deformations of the LF and the DNLS 
equations are proposed. 
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I. INTRODUCTION 
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Intensive research extended over fifty years, starting from the pioneering work of GGKM ggkm65 
has taken the theory of nonlinear integrable systems to a well established subject, seemingly with 
no further surprises could be expected from it. However, quite recently an integrable 6th-order Ko- 
rteweg de Vris (KdV) equation was discovered 6kdv, apparently contradicting the conventional wisdom. 
Though eventually this was understood as an nonholonomic deformation of the standard KdV equa- 
tion kuper08,kundu082,chineese08, this unconventional concept stimulated further study to unveil the 
nonholonomic deformation for the entire AKNS family JMP09. Consequently, deformed integrable hier- 
archies were discovered for the nonlinear Schrodinger (NLS) equation, the sine-Gordon equation and the 
modified KdV equation JMP09,kundu082, kundu083, though the basic idea, exploiting the negative flow 
in the spectral parameter, was already used in some specific cases stramp,dmkdv09. Such deformations 
for the field theoretical models can be interpreted also as a perturbation of the original integrable system 
with certain differential constraint on the perturbing function, such that the perturbed system together 
with their deformed hierarchy as a whole becomes a new integrable system. However, in spite of the nat- 
ural expectation JMP09, the nonholonomic deformation could not be constructed for the Kaup-Newell 
(KN) class KN, which includes the important integrable system like the derivative NLS (DNLS) equation 
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In such a situation, another unexpected integrable equation (|2ip was proposed very recently by Lenells 
and Fokas (LF) LF, soliton solution of which exhibits close resemblance with the solution of the DNLS 
equation. This intriguing fact gave us the right hint to formulate the nonholonomic deformation for 
the KN spectral problem and applying this construction discover a novel mixed integrable hierarchy for 
the deformed DNLS family, yielding more general accelerating soliton solutions. This unusual soliton 
dynamics is sustained by the time-dependent asymptotic of the perturbing functions, which themselves 
take a consistent soli tonic form. Solving the differential constraints it is possible in some cases to express 
the perturbing functions locally through the basic field, reducing the perturbed system to a novel higher 
nonlinear equation in terms of the basic or its potential field. 

Most interestingly, we could identify the LF equation as the first nontrivial member in the mixed 
hierarchy of the deformed DNLS, which enables us to bypass the lengthy and involved approaches like 
the Riemann-Hilbert (RH) or the dressing method for extracting the soliton solution of the LF equation 
adopted in LF,lenells09 and obtain the same result by suitably deforming the well known DNLS soliton. 
Moreover, such LF solitons can also exhibit an unusual accelerating motion due to the time-dependent 
boundary condition of the deforming functions. 

Considering the deformed DNLS hierarchy we can also construct a novel higher order deformation of 
the LF as well as the DNLS equation. 

As a logical extension we discover the nonholonomic deformation for another form of the DNLS 
equation (06]), proposed by Chen-Lie-Liu (CLL) c]l, which is found to share the similar soliton dynamics 
and the mixed deformed integrable hierarchy, mentioned above. 

The arrangement of the paper is as follows. Sec. II presents the known result of the undeformed KN 
system. Sec. Ill formulates the nonholonomic deformation for the KN family, including its deformed 
two-fold integrable hierarchy. Sec. IV classifies these deformed hierarchies, identifying the LF and the 
deformed DNLS equations. Sec. V presents the details on the LF equation as a member of the deformed 
KN hierarchy together with its higher deformations. Sec. VI details similar result for the deformed 
DNLS equation. Sec. VII derives in a simple way the exact accelerating soliton solution of the deformed 
equations, including that of the LF equation. Sec. VIII introduces the CLL equation, its nonholonomic 
deformation and the related soliton solution with acceleration and other peculiarities. Sec. IX is the 
concluding section, followed by the bibliography. 



II. KN SPECTRAL PROBLEM FOR THE DNLS EQUATION 

For solving analytically the DNLS equation 

iqt ~ Qxx + 2ie(\q\ 2 q) x = 0, (1) 

through the inverse scattering method (ISM), Kaup and Newell introduced a new type of spectral problem 
given by the Lax pair KN: 

U dnls (X) = -i\ 2 a 3 + XU^, = qa + + eq*a-, (2) 

WA) = (ia 3 U^ + (U^f)X-ia 3 (U^) 2 X 2 + 2U^X 3 -2ia 3 \ 4 . (3) 

The flatness condition 

U t -V x + [U,V] = (4) 

of this Lax pair (|2l3p generates the DNLS equation (P) together with its complex conjugate. Note that 
the time-Lax operator ([3]), associated with the second order dispersive DNLS equation, contains positive 
powers up to N+ = 4 in the spectral parameter A. The higher order integrable equations in the DNLS 
hierarchy can be obtained similarly from a Lax pair like (0]) , where the space-Lax operator U(X) remains 
the same, while the time-Lax operator V(A) is replaced by V^+)(A), containing higher positive powers 
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in A. For example, for generating a Z-th order dispersive nonlinear equation in the DNLS hierarchy, the 
V^ N+ '(X) matrix should contain powers in A up to N+ = 21. 

Recall that for extracting the solutions of the DNLS equation through the ISM KN, the key role 
is played by the space-Lax operator ([2]), the explicit knowledge of which together with the required 
vanishing boundary condition for the fields are sufficient for constructing the exact form of the soliton 
solution. Only at the final stage of the solution we need to know the asymptotic value of the time 
evolution operator ([3]): Vdnis{X)\\ x \^oo = ^oo(A) = — 2i\ 4 a 3 , for determining the time-dependence of the 
spectral data S = {p(X 2 ), for A real;^, for A^ discrete, k G [l,iV]}. Pure soliton solution is obtained 
by setting the reflection coefficient p(X 2 ) = 0, when the time evolution of the spectral data is obtained 
as Afc = 0, 6^ = — 4i\fcbk(t) yielding Xk(t) = A/ c (0), bk(t) = e _4 * A fc*6fc(0), which in turn fix the soliton 
velocity vq and the frequency ojq of the wave enveloping the soliton. 

Note that for all higher order equations in the DNLS hierarchy the entire picture of the ISM, 
including the structure of soliton solutions basically remains the same. Only the nontrivial time- 
dependence of the spectral data, determined by the boundary condition of the associated time-Lax 

operator Vj^^ = —2iX N+ a i becomes more general: bk(t) = e~ ilX * +t bk(0), changing the soliton param- 
eters vo,u>o, accordingly. 

Referring to KN for details on the N-soliton solution, we reproduce here the well known 1-soliton of 
the DNLS equation (P) as 

a mn o- n exp(46>) + exp(-ie7) 

q = 4er? exp(26> - 2ia) . /A ' ^ — ^ (5) 

v y (exp(4(9) +exp(ie7)) 2 w 

where 

a = ix + 2uj t, e = 7](x- 2v t), X\ = £ + if] = AV 7 (6) 

with ojq as the frequency of the enveloping wave and vq as the constant soliton velocity. We shall see 
that the soliton solution of the nonholonomic deformation of the DNLS hierarchy of equations can be 
obtained in a easy way by deforming suitably the undeformed soliton ([5]). 



III. NONHOLONOMIC DEFORMATION OF THE KN SYSTEM AND ITS MIXED IN- 
TEGRABLE HIERARCHY 

Recall that, perturbation as a rule destroys the integrability of a system. The idea of perturbation 
through nonholonomic deformation however is to perturb an integrable system with a deforming function, 
such that under suitable differential constraints on the perturbing functions the integrability of the 
whole system is preserved again. In a field theoretic models like the present one, a constraint given 
by differential relations (not evolution equations) is equivalent to a nonholonomic constraint. For such 
novel nonlinear systems, like the perturbed DNLS equation, we can derive all the important features of 
a genuine integrable system, i.e. the Lax pair, the integrable hierarchy, the exact N-soliton solution etc. 

For implementing the concept of a nonholonomically deformed system we start not from the perturbed 
equation, but by constructing the associated deformed Lax pair U (A), V(X), the flatness condition of which 
generates then the required integrable perturbed equation. The space-Lax operator for this perturbed 
system is chosen to be the same as that of the undeformed operator Uq(X), while in the time-Lax operator 
a deforming part should be added: V(X) = Vo(X)+Vd e f(X). Therefore, since we keep the spectral problem, 
central to the ISM, unchanged the deformed systems automatically become exactly solvable by the ISM. 
However due to the deformation of the time-Lax operator, the time evolution of the spectral data gets 
changed. As a result one gets a deformed soliton solution with changed dynamics, where interestingly, 
along with the basic field the perturbing functions also take similar solitonic form. A point of significant 
practical importance in such perturbed integrable equations is that, though the ISM or other exact 
methods like the RH method, the dressing method, Hirota's bilinearization etc. are applicable here, one 
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can extract in fact the same result bypassing these complicated methods. For this one has to take the 
well known soliton solution of the unperturbed field and change its dynamics by simply deforming the 
soliton velocity and the frequency of the enveloping wave, which can induce more complicated soliton 
dynamics, including accelerated or decelerated soliton motion, as we will see below. 

A. Nonholonomic deformation of the DNLS and the perturbed integrable hierarchy 

For generating the perturbed integrable systems through nonholonomic deformation of the DNLS 
equation, we take the space-Lax operator same as in the unperturbed case ([2]), while for the time-Lax 
operator there could be different choices. We may start from any known integrable hierarchy of the DNLS 
associated with Vo(A) = V^+^A), which represents a polynomial in the spectral parameter X n with its 
power n = 0, 1, . . . running upto an arbitrary positive even integer N+ . Such higher operators can be 
constructed systematically in the form 

N+ a i 

V( N +\x)=Y,W k (^-) l (U {0) ) m , = qa + + eq*a-, N + = -(k + m) + l, (7) 
dx 2 

where the summation indices k,l,m > are all positive integers with the total 1 + iV? number of 
terms appearing in this expression JMP09. Therefore the corresponding unperturbed system represents 
an integrable higher order equation in the +ve hierarchy of the DNLS equation with iiV_|_-th order 
dispersion. Since we are considering the deformed system, the related time-Lax operator is constructed 
as V(A) = Vo(A) + Vdef(X), as stated above, with different options for the deformed Vd e f(X) operator 
containing negative powers of A : 

N— 

v£f\x) = iJ2x-iGW. (8) 

j=0 

G^\ j = 0, 1, 2, . . . iV_ is a set of deforming matrices with arbitrary even integer iV_ , higher values 
of which indicate higher orders of deformation. It is important to note that, while the Lax operators 
U(X) and the undeformed part Vq(X) involve only the basic field q, q\ the deforming part vj^i (X) 
contains only the perturbing functions, expressed through the matrix elements of the deforming matrices 
j = 0, 1, 2, . . . iV_. The nonholonomic constraints on such perturbing functions are represented by 
coupled differential equations, defined in turn by the flatness condition of the Lax pair for the deformed 
system. With the above arrangement we can build the deformed KN system again through the condition 
(jll) as 

Uot - V 0x + [Uo, Vo] = V defx - [U , V def ], (9) 

yielding at different powers of the parameter A the perturbed equation together with the nonholonomic 
constraints on the deforming matrix Vdef ■ Note that though the basic field equation obtained from ([9]) may 
differ as different members of the positive flow hierarchy depending on the choice of ./V+ in Vo = V^ N +\X), 
the perturbative part with deforming matrices , j = 0, 1, 2, . . . , N_ with the nonholonomic constraints 
would independently yield another hierarchy of recursive relations 

G (o) = pW^-ifaG®], 

G (D = [u«»,GW}-i[a 3 ,G(% ... , 

= i[uV\G( N -V]-i[a 3 ,G( N -\ 

Gf-D = i[tf(°),G*"->], Gf-)=0, (10) 

obtainable also from the flatness condition of the deformed Lax pair. Remarkably, this deformed KN 
hierarchy differs considerably from the deformed AKNS hierarchy found in JMP09 due to different spectral 
dependence between the corresponding Lax operators U(X) and the appearance of an additional deforming 



4 



matrix in the KN case, which leads also to a peculiar G^ N ~\t) matrix depending only on the time 
variable t, as seen from the last equation in (jlOl) , 

Using the explicit form of the Lax operators from ([2]) and ([3]) or its higher hierarchical form (|7|) 
together with the deforming operator ([8]), one can generate explicitly a two- fold integrable hierarchy of 
the perturbed KN-DNLS family, which we present below in detail. 

IV. CLASSIFICATION OF THE INTEGRABLE DEFORMATIONS IN THE KN HIER- 
ARCHY 

Before concentrating on the soliton solutions and detailed structure of the deformed DNLS and the 
LF equations in the subsequent sections, we present here a classification for the deformed KN hierarchy. 
For this we take the space-Lax operator in the standard form ([2]), but select a more general time-Lax 
operator 

V^+' N -\X) = V^ N+ \\) + v£f\\) (11) 

with mixed flows in the spectral parameter exhibiting a two-fold integrable hierarchy. This includes a 
positive flow with 

y( JV +)(A) = V + XV l + --- + X N +V N+ , (12) 

representing the known higher order integrable hierarchy of the DNLS equation, together with a negative 
flow: 

v£f\\) = i(G<°> + A- 1 ^ 1 ) + • • • + \- N -G^), (13) 
describing the deformation hierarchy of the DNLS equation given through perturbation matrices ([8]). 

A. Hierarchy without deformation 

This is generated by the pure positive flow with arbitrary even integer N+ and N- = 0, having the 
time-Lax operator: y(-^+'°) = y( N +) i as given in (|12p . This case corresponds to the well known KN- 
DNLS integrable hierarchy without deformation KN, where each of the constituent matrices Vj in the 
expansion 

y(N + ) = X jy. 

can be build up systematically as ([7]) from the known building blocks using 
only a dimensional argument, as shown in JMP09. The particular choice N + = 4 yields the standard 
undeformed DNLS equation. 

B. Pure deformation hierarchy 

This is the complimentary case generated by the pure negative flow with N + = and iV_ as an 
arbitrary even integer with the time-Lax operator: V^ 0,N ~^ = vj^T , given solely through deformation 
matrices (I13D . Such deformations can be interpreted also as novel integrable perturbations of the lowest 
order equation qt = in the KN hierarchy, where we can generate perturbed equations with different 
nonholonomic constraints for different choices of iV_ in (|8|). We find remarkably that, any deformation 
with iV_ < 2 leads to a trivial result with the lowest order deformation produced at iV_ = 2 and all 
the subsequent higher order deformations seem to be given by higher even numbers iV_ > 2. We would 
analyze below in detail such deformations with concrete examples. 

C. General deformed case with mixed flow exhibiting a two- fold hierarchy 

This mixed integrable hierarchy is generated by both positive as well as negative flows with N+ 
and iV_ being arbitrary even integers, where the time-Lax operator is given in the general form (1 1 1 1> : 
y(N+,N-) = y(N+) + V^f\ containing CGD as well as (fl3"|) . The known integrable KN hierarchy of KN is 
perturbed through another novel deformed hierarchy, producing thus a two-fold hierarchy of the deformed 
DNLS equations for different choices of N+, iV_, few interesting cases of which we present below. 

i). Lenells-Fokas equation: 
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The first member with both positive and negative flow of the KN system is given by the lowest 
nontrivial values N+ = 2 and iV_ = 2, with the time-Lax operator: y( 2 > 2 ) = + vj^J. Note that from 

the dimensional argument JMP09 we must have ~ U(X) = all(X) and vj^J should be constructed 
from (|13p with minimum number of deforming matrices G^\j = 0, 1, 2. Surprisingly, one finds that the 
recently proposed LF equation LF is given in fact by this mixed deformed KN hierarchy with the lowest 
order deformation, where N- > 2 would generate a novel integrable hierarchy of the LF equation with 
higher order deformations. We present the case of LF equation in detail in the next section. 

ii). Deformed DNLS equation 

The next entry in the mixed deformed KN hierarchy is given by N + = 4, iV_ = 2, with the time-Lax 
operator: l/( 4 ' 2 ) = vj^ s + V^j, with vj^ ls as the well known DNLS operator ([3]) and the deforming 

operator V^L same as in the above Lenells-Fokas case constructed from (|13p with N_ = 2. Therefore, 
since iV_|_ = 4, iV_ = corresponds to the unperturbed DNLS equation, the present case with N+ = 
4, N- = 2 would produce its integrable perturbation with minimal nonholonomic constraint on the 
deforming functions. Higher even integer values of iV_ > 2 would generate the DNLS hierarchy with 
higher order integrable deformations. 

We study below different cases of the deformed KN hierarchy in detail. 

V. LENELLS-FOKAS EQUATION AND ITS DEFORMED HIERARCHY 

Recently a new type of integrable NLS equation is proposed by Lenells and Fokas and solved it exactly 
for the soliton solution through the RH method LF. They also observed certain similarity of this solution 
with the well known DNLS soliton KN. We find intriguingly that, the LF equation can be given by 
the minimal nonholonomic deformation of the first member in the KN hierarchy, where the subsequent 
perturbation can generate a novel LF hierarchy with the higher order deformation. The link of the LF 
equation with the first conservation law in the negative flow of the KN hierarchy has also been identified 
quite recently by Lenells lenells090. These findings resolve thus the mystery around the integr ability 
of the LF equation and the resemblance of its solution with the DNLS soliton. As a bonus we can get 
also the explicit soliton solutions of the LF equation in an easy way by deforming the well known DNLS 
soliton and thus can derive the same result of LF,lenells09 through much simpler route, bypassing the 
involved RH or dressing method employed in LF,lenells09. Moreover the asymptotic form of our time-Lax 
operator given through the time-dependent boundary condition of the deforming functions can lead to a 
remarkable possibility of creating accelerating solitons of the LF equation, similar to that found for the 
deformed KdV equation kund j u082. Following the scheme described above we can construct along with 
the LF equation a novel integrable LF hierarchy with increasingly higher order deformation. 

The associated Lax pair for the LF equation, as classified above for the deformed KN-DNLS hierarchy, 
can be constructed from the space-Lax operator U(\) of the DNLS equation ([2]) and the time-Lax operator 
given through deformation of V^ 2 \X) = all(X) as 

V^ 2 \X) = aU(X) + vQ(X), 
V £}(X) = i(G<°> + A- 1 ^ 1 ) + X- 2 G^) (14) 

It is intriguing that with a variable change (x, t) — > (x,t — ^x) one can absorb the positive flow part 
V"( 2 )(A) = aU(X) from the time-Lax operator (|14p . simplifying the LF equation and bringing this system 
to a pure deformation case discussed in Sec. IV B. We can derive the constraint equations on the deforming 
matrices G^\j = 0, 1,2, from the general form (|10p for this minimal deformation with iV_ = 2 as 

G^ = [U^,G^]-i[as,G(-% GW = [U(°\G(% G® = 0, (15) 
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where = qa + eq*a + as in ([2]). Consistent with these relations we may express the deforming 
matrices as 

G(°) = g( x , t)a 3 , = i(w(x, t)a + + ew(x, t)*a~), G {2) = a(t)a 3 , (16) 

through the deforming functions g(x, t),w(x, t) and an arbitrary function a(t), dependent only on time t. 
We can derive the perturbed field equation from the relation ([9]) by using the deformed Lax pair (pi I14p : 

[/f> - a uj® = -i[U^\G^\+i[a\G^\. (17) 

Here the original equation is the LHS part, while the perturbation through deforming matrices is given 
by the RHS terms. Using now the structure of J/(°) and G^,G^ as in (fT6|) . we obtain the perturbed 
field equation from the matrix elements of (|17p as 

q t - aq x + 2w + 2igq = (18) 

together with its complex conjugate. Note that here the nonlinear equation is generated by perturbation 
only, from an original linear equation qt — aq x = 0. In this scheme w and g are like perturbing functions 
which are subjected to the nonholonomic constraints obtained from (|15p as 

g x (x,t) = ie(qw* — q*w), w x (x, t) = — 2ia(t)q. (19) 

Note that we can consider equations (|18p as a perturbed integrable equation with the nonholonomic 
constraints (|19p imposed on the perturbing functions. However since the original equation in this case is 
only a nonsignificant linear equation, an interesting alternative approach here would be to obtain a new 
nonlinear integrable equation for the basic field by resolving the constraint relations and expressing all 
perturbating functions through this field itself. Therefore for simplifying the set of equations (|18til9[) we 
introduce a potential field q = u x , which solves the constraints on the perturbing functions through the 
potential field as 

g(x,t) = -2ea(t)\u\ 2 + c(t), and w(x, t) = -2ia(t)u, (20) 
and leads (|18p finally to an evolution equation with a nonlinear derivative term in the form 

u x t — au xx + 4ia(t)u + 2ic(t)u x — 4iea(t)\u\ 2 u x = (21) 

together with its complex conjugate and leaving no more constraint. We note immediately that (|2ip is a 
nonautonomous, though integrable generalization of the LF equation LF with arbitrary time-dependent 
coefficients a(t),c(t). For constant choices of these functions: a(t) = ao,c(t) = cq, on the other hand 
we get the corresponding autonomous LF equation and for a particular choice ao = l,co = we recover 
immediately from (|2ip the recently proposed LF equation 

u x t — au xx + 4iu + 4ie\u\ 2 u x = 0. (22) 

It should be further noted that, a simple change in time- variable: t — > t — ^x, as mentioned above, can 
remove the dispersive term —au xx from the LF equation (|22p bringing it to a simpler form (|24p considered 
below. We show in Sec. VII. that for the time-dependent a(t) = aot and c(t) = cot the more general LF 
equation (f2Tj) allows novel accelerating soliton solution. 

A simpler LF equation 

A simpler equation equivalent to the LF equation, representing the DNLS hierarchy with pure defor- 
mation can be derived from the general case at iV + = 0, iV_ = 2. The Associated Lax pair is given by a 
similar form (|14p . with the term V^ 2 '(X) removed, i.e along with ([2]) one should have V^ 2 \\) = V^j(X) , 

where the deformed Lax operator V^ e i(A), is the same (114p as given for the LF case. We have mentioned 



7 



already, that the LF equation can be brought to this simpler case with the time-Lax operator V^°' 2 \\) 
by a simple change in the time- variable. As a result we get the set of equations given by (|18h without 
the —aq x term: 

q t + 2w + 2igq = 0, (23) 

together with its complex conjugate. The constraint relations (I19D however remain the same, leading to 
a unconstraint nonlinear integrable equation 

u xt - 4i(e\u\ 2 u x - u) = 0, (24) 

where we have taken a(t) = 1, c(i) = as for the LF equation. This simpler integrable equation (I24j) . 
which is equivalent to the LF equation under a change of the time variable, is the least possible pure 
deformation of the DNLS hierarchy. Soliton solutions of this simplest member of the deformed KN-DNLS 
equations can be found similarly by deforming the known DNLS soliton. 

B. Novel higher order deformation hierarchy of the LF equation 

For generating the deformed hierarchy of the LF equation given by the higher order nonholonomic 
constraint on the perturbing functions, one has to take iV_ as an higher even integer, while keeping N + = 2 
fixed as in the LF equation. This would lead to the same perturbed equation (fT7|) or equivalently to (fTBj) . 
though due to the increase in the number of deforming matrices from j = 0, 1, 2 to j = 0, 1, 2, . . . , A_ 
the constraint equations (|19p would be replaced by the more general deformed hierarchy (|10j) . Therefore 
it is no longer possible to solve the perturbing functions through the potential field in a easy way like 
(I20p . However it should be a challenge to derive some generalization of the LF equation (|2ip for this 
higher deformed hierarchy expressed again through the potential field u. Therefore we take up the case 
of the first higher deformation with iV_ = 4, assuming e = 1 for simplicity. Here a larger set of deforming 
matrices j = 0, 1, 2, 3, 4, can be expressed as 

= g( x , t)a 3 , = i(w(x, t)a + + w*(x, t)a~), G (2) = a(x, t)a 3 
G (3) = i(b(x,t)a + + b*(x,t)a~), G (4) = c(t)a 3 (25) 

through the deforming functions g(x,t),w(x,t),a(x,t),b(x,t) and an arbitrary time-dependent function 
c(t). Using the nonholonomic constraints (|10p for iV_ = 4 and introducing the potential field q = u x 
as above, we can partially remove the constraints and get a novel higher order deformation of the LF 
equation (i22j) as 

u x t — olu xx + Aiu + Aie\u\ 2 u x + 2igu x + 2w = 0, (26) 
where the constraints on the residual perturbing functions g, w are given by 

9x = i(u x w* - u* x w), w x = 4c(t)(i\u\ 2 u x - u). (27) 

Note that the first part of the deformed LF equation (|26p is exactly same as the original LF (|22p . which 
however is deformed now by the effective perturbing functions, with new nonholonomic constraints (|27p . 

We can construct the soliton solutions of the LF equation as well as all its higher deformations using 
a simple trick of deforming the well known DNLS soliton, which we demonstrate in a later section. 



VI. NONHOLONOMIC DEFORMATION OF THE DNLS EQUATION AND ITS DEFORMED HIER- 
ARCHY 

We concentrate now on a more interesting deformed DNLS equation, which may be obtained from 
the general Lax pair (jlip for the particular choice = 4, A r _ = 2. Therefore in this case the Lax 
operator U(X) is the same as ([2]), while V(X) as mentioned already is given by 

^ (4 ' 2) (A) = ^i 4 L(A) + ui 2 ;(A) (28) 
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where V^ nls (X) is the well known time-Lax operator ([3]) of the DNLS equation and V^ e j(A) is the deforming 
part (I14p . same as in the LF case. The flatness condition of this Lax pair yields the deformed equation 
in the matrix form; 

Of) _ [/(O) _ ((U^f) x = i(-p7< >,G<°>] + [a 3 ,G^]), (29) 

where the matrix U^ ' is defined through the basic fields q, q* as in ([2]) and the nonholonomic constraints 
on the perturbing matrices are given again by (|15p . In the scalar form the perturbed DNLS equation 
(|29l) may be expressed as 

iqt + Qxx ~ 2ie(\q\ 2 q) x = 2(gq - iw), (30) 

where the LHS is the original DNLS equation, while the RHS gives its integrable perturbation through 
functions w,g , with the nonholonomic constraints on them defined as in (|19|) . 

We may resolve these constraints by introducing the potential field q = u x , and express the perturbing 
functions through potential field as (|20p . same as in the LF case. Inserting this set of relations in the 
deformed DNLS (|30|) . we obtain a new integrable potential DNLS equation in the form 

iu xt + u xxx - 2ie(\u x \ 2 u x ) x - 4a(t)u - 2c(t)u x + 4ea(t)\u\ 2 u x = 0. (31) 

This equation generalizes the LF equation (|2ip and its simpler version (|3ip with the addition of another 
nonlinear derivative term together with a higher dispersive term. Notice that considering t — > it and 
u = iuq with «o as a real field, one can obtain from (|31|) an integrable equation, which may be interpreted 
as a novel derivative generalization of the mKdV equation. 

Deformed hierarchy of the DNLS equation may be constructed following the scheme for the LF system 
considered above, with the only difference that, in the DNLS case one has to take N + = 4 as in ([3]), 
while iV_ for both these systems may be any even integer as in (|13p . This results again to the same 
perturbed DNLS equation (|30p . though the nonholonomic constraints should be given now by a more 
general deformed hierarchy (I10p . For iV_ = 4, which corresponds to the first higher deformation similar 
to that considered above for the deformed LF equation, we obtain the same perturbed DNLS equation 
(|30p . though with higher order nonholonomic constraints on the perturbing functions as 

g x = i(qw* - q*w), w x = -2i(aq + b), a x = i(qb* - q*b), b x = -2ic(t)q, (32) 

where the functions g, w, a, b and c(t) are the elements of the five deforming matrices (|25p . 



VII. EXACT ACCELERATING SOLITON SOLUTIONS OF THE DEFORMED DNLS 
HIERARCHY 

As we have seen through the examples of the deformed DNLS including the LF equation, for con- 
structing the nonholonomic deformations one has to take the deformed time-Lax operator V(A) with 
an additional part Vd e f(X), while keeping the same space-Lax operator U(X) JMP09. As a result the 
x-dependent part in the deformed soliton solution remains the same as in the undeformed case, while its 
time-dependent part gets changed by bringing in a new time in its dynamics, defined by the asymptotic 
value of lim x ^± 00 V(A). Interestingly, since the ISM for finding the soliton solutions is based primarily on 
the spectral problem $ x = U(A)<& associated with the space-Lax operator, the structure of the N-soliton 
solution remains the same even for all the deformed DNLS equations in their two-fold hierarchy with 
arbitrary N + , N- . Only at the final stage of the ISM we need to express the time-dependent part of the N- 
soliton by inserting the time-dependence of the related discrete spectral data: Xj, bj(t), j = 1, 2, . . . , N, 
determined through 

N- N- 

^lirn^ V(X) = i(2X N + + A _2i c 2 Kt)) v 3 , as A fc = 0, b k (t) = i{2X% + + £ X^ l c 2l (t)) b k (t). (33) 
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Here the time-dependent arbitrary functions C2i(t), I = 0, 1, 2, • • • are the boundary conditions of the 
deforming functions \\m. x ^± 00 G^ = icj(t)a 3 , j = 0, 2,4, . . . iV_. Note that, in (|33|) we have restricted 
only to the terms with even powers in A, since it is evident from the structure of the deforming matrices 
(|16p and (|25p in the two cases of deformation with N- = 2 and N_ = 4 that we considered above, that 
only the matrices G 3 with even indices j = 0,2,4 are diagonal and may have nontrivial time-dependent 
asymptotic at the space-boundaries, while those with odd indices vanish at space infinities with their 
corresponding perturbing functions along with the basic field. Note that the first term in (133p comes 
from the standard hierarchy with positive flow V^+^A), while the rest of the terms is due the negative 
flow in this two- fold hierarchy, which describes the evolution along a new time, generated by (133 j) . 

Therefore we can simplify the solution procedure enormously, since instead of applying the ISM or any 
other sophisticated method for extracting the soliton solution individually for every equation belonging 
to this two-fold infinite hierarchy, we can simply take already available well known N-soliton solution 
of the original DNLS equation as in Kaup-Newell KN and deform the time-dependent part by suitably 
choosing the values of N± and thus obtain the N-soliton solution for any deformed equation in this 
integrable hierarchy almost without any effort. For example, using the identification of the Lenells-Fokas 
equation LF as the N + = 2, iV_ = 2 member in the deformed DNLS hierarchy, we can reproduce its 
soliton solution quite easily, without going trough the involved and lengthy procedure of the RH or the 
dressing method used in LF,lenells09. Moreover choosing the time-dependent boundary condition of the 
deforming fields suitably we can find also a novel accelerating soliton in such deformed KN systems, 
observed already for the deformed AKNS JMP09. 

Based on our above argument we present the well known 1-soliton solution ([5]) of the KN-DNLS 
equation by deforming its parameters u — * w(t), v — * v(t) KN 

where 

a = ix + 2w(t), 8 = rj(x- 2v(t)), \\ = £ + i-q = AV 7 (35) 

Such deformations can be extended easily to the N-soliton solution, which however we will not detail 
here. We stress that the 1-soliton solutions of all the equations in the deformed DNLS hierarchy can be 
expressed by the same expression (|34|) . only the form of the enveloping frequency oj(t) and the soliton 
velocity v(t) would be different for particular equations, depending on the time-dependence of the 
corresponding spectral data, as obtained from the general situation ([33]) . 

For example in case of the standard DNLS KN with N + = 4, A r _ = 0, the time dependence obtained 
from ([33]) is lim x ^± oc Vdnis{^) = — 2iA 4 cr 3 , which gives therefore u>(t) = oj^t and v(t) = vot, where constant 
parameters loo,vq are given in as found through the ISM in the pioneering work of Kaup-Newell KN. A 
similar though more generalized formulation, giving the possibility of accelerating solitonic motion, will 
be presented below for the nonholonomic deformations of the DNLS hierarchy including the LF equation. 

A. Accelerating soliton in the LF equation as deformation of the DNLS soliton 

For the LF equation (|31~1) the time dependence of the spectral data, as explained above, is to be 



determined from ([33} with N + = 2, JV_ = 2 as bi(t) = exp{-2i(a\jt - J dt(c(t) + ^B)))h(0). Therefore 
following the above argument we can derive the soliton solution for the LF equation from (|34p quite easily 
by specifying the time-dependence of the parameters with deformation as 

rt a u\ 

u(t) = 2ite($(Ai)) = -2atA 2 cos7 + 2 / dt'(-^f cos 7 - c(t)), 

2 r f c(t') a(t') 

v(t) = -Im^(X 1 )) = -2at-2 j dt'( -U- sin 7 - -i^), 
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where $(A a ) = a\\t - f dt'(c(t') + ^Q-). 



(36) 



Here the arbitrary time-dependent functions c(t), a(t), which controls the dynamics of the soliton from the 
space-boundaries are given by the asymptotic values of the deforming matrices and , respectively. 
Thus we can obtain the soliton solution (|34p for the LF equation with the time-dependent soliton velocity 
v(t) and frequency oj(t) as (|36p . which can yield accelerating solitons, generalizing the solution found in 
LF,lenells09. In particular, for the choice of the functions c(t) = cot, a(t) = dot we obtain an accelerating 
LF soliton: 

a (no o- x exp(46>) + exp(-ie7) 

4e7/exp(2# — 2ia)- — — - — --^ (37) 



(exp(46>) + exp(ie7)) s 



where 



d = £,X - LU t + LOdt , 
= T](X ~ (V t - V d t 2 )), 



where ujq = 2a A 2 , ujd 



-((^) cos 7 + co), 



where v = 2a, v d = sin 7 - 



(38) 



Integrating the solution (I37p we can get the soliton solution of the LF equation (|2ip for the potential 
field moving with a constant acceleration as 



-2ia(t) sin 7 



exp(— (icq + 2ia)) 
exp(26») + exp(ie7- 26)' 



(39) 



It is evident that for Vd(t) < we get an accelerating, while for Vd(t) > a decelerating soliton. Such 
a soliton for the generalized LF equation (|2ip for the modulus of the potential field 



\u(x,t)\ 



2 sin 7 



(cosh 2 29 - sin 2 2 



(40) 



given through solution (|39p is presented in Fig. 1. Note that (I39p gives a more general soliton solution 




Figure 1: Dynamics of the exact soliton solution (|39p for the potential field \u(x, t)\ of the general- 
ized Lenells-Fokas equation ([2T]) . The unusual shape change and bending of the soliton indicates its 
decelerating motion with time-dependent velocity Vd(t). 

of the LF equation, where for the particular choice of the parameters with constant values: c(t) = 
and a(t) = 1, one recovers exactly the constant velocity soliton found in LF. This demonstrates that we 
can extract the same soliton solution of the LF equation by just deforming the parameters of the known 
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soliton KN of the DNLS equation and thus can avoid the use of sophisticated methods like the RH or 
the dressing method LF,lenells09. 

Another significant achievement is to construct in the same simple way the soliton solutions for all 
higher deformed LF equations by an easy extension of the time-dependence of the known DNLS or the 
LF soliton. For example, for the first higher deformation of the LF equation given by the higher nonlinear 
extension (|26l) . the soliton solution can be expressed again in the form (139 j) . with similar expressions for 



a = £x + 2u(t), 9 = rj(x- 2v(t)) 
where the time-dependence of the parameters should be given now in a more extended form 

uj{t) = 2ite($(Ai)), v(t) = -/m($(Ai)), 

V 

where $(Ai) = aXjt - f dt'(g (t>) + ^ + ^), (41) 

in comparison with the LF soliton (|36p . The arbitrary time dependent functions go (t) and ao(t) are 
asymptotic values of the perturbing functions g(x, t) and a(x, t) as in ([25]) . respectively and can naturally 
be taken also as constant parameters. 

Similarly, the N-soliton solution of the LF equation can be obtained avoiding the involved methods 
of LF,lenells09, by deforming the known N-soliton of the undeformed DNLS equation KN using (|33|) 
with the choice N + = 2,N— = 2, which corresponds to the LF equation. We are however not going to 
reproduce here these otherwise simple steps. 

B. Accelerating soliton in the deformed DNLS equation 

As we have stated, the soliton solutions for the integrable nonholonomic deformation of the DNLS 
equation can be derived using general techniques like 1ST, Hirota's bilinearization etc. or by the RH 
and the dressing method as done in LF,lenells09. However we find them here in the simplest way by 
deforming the well known DNLS solitons already obtained in KN through the ISM. We start therefore 
from the known undeformed soliton (|43p and deform just its parameters like frequency and the velocity 
using (|36p or (138 j) . as we done for extracting the LF soliton. Thus we derive the 1-soliton for the deformed 
DNLS equation (|30til9|) or equivalently for (|3ip . again in the form (|37p yielding 



(cos 8 + cos 7) 



with the only difference that the undeformed parameters like the initial frequency and the velocity 
are given now as in the original DNLS soliton KN: 

u = 2Re{\\) = 2{f -if) = 2A 4 cos 2 7 , v = 2Im{\\) /rj = 4£ = 4A 2 cos 7. (43) 



The deformed DNLS soliton (|42p as shown in Fig. 2a may exhibit again accelerated or decelerated 
motion, similar to the LF case depicted in Fig. 1, but with a difference in their initial velocities Vq. It is 
important to note here, that the perturbing functions g, w, which maintain the integrability of the system 
are not any given functions as in the usual perturbation theory, but are determined consistently through 
internal dynamics and can be expressed for the present models through the potential field u(x, t) as (|20p . 
Therefore for the soliton solution (|39p of the basic field these perturbing functions also take the solitonic 
form, which for the LF and the deformed DNLS equations in particular is as 

w(x,t) = 2ia(t)u (44) 
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expressed through the solution ([39]) for the potential field u(x,t) and similarly using ([40]) : 



g(x, t) = 2ea(t)\u\ 2 + c(t) = c(t) + 8ea(t) 




(45) 



cosh 2 29 - sin 2 \ ' 



graphical representation of which is given in Fig. 2b. Note that w(x,t) is a complex function with soliton 
solution (|44p . while g(x,t) is a real function with solution (|45p . both having time-dependent velocity 
v(t) and frequency co(t) as in ([38]) . but additionally with a time-dependent amplitude a(i). It should be 
noticed that, while at space infinities the perturbing function w vanishes, g goes to a time-dependent 
asymptotic c(t), and drives the solitons from the space-boundaries. In this way the perturbing functions, 
themselves being consistent soliton solutions, preserve the integrability of the system and can force the 
basic field soliton to accelerate or decelerate. Such perturbations thus can play potentially important 
role in practical applications as we discuss later. 



Figure 2: Exact accelerating solitons in the nonholonomic deformation of the DNLS equation, a) Soliton 
for the basic field \q\ corresponding to (|37p with (|43p . b) Consistent solitonic form (|45p for the perturbing 
field g(x,t) expressed through the potential field \u(x,t)\ 2 . Notice the nontrivial boundary value of the 
perturbed function. 

The higher N-soliton of the nonholonomic deformation of the DNLS equation can similarly be derived 
by deforming the known N-soliton of the standard DNLS equation KN. Likewise the solutions of all 
integrable equations in the deformed hierarchy can be obtained from those of the original undeformed 
integrable hierarchy by suitable deformation of the solitonic parameters as described above. In such 
deformed hierarchies, for every fixed N + , one may choose N_ as any even integer values, like iV_ = 2,4, 
as we have considered above, generating thus a novel two-fold hierarchy. 

VIII. NONHOLONOMIC DEFORMATION OF THE CLL EQUATIONS, ITS TWO-FOLD 
HIERARCHY AND ACCELERATING SOLITONS 

The CLL equation is another form of integrable derivative NLS equation with a different type of 
nonlinearity ell. Following our approach as applied above for the deformed DNLS equation, we construct 
here the integrable nonholonomic deformation for the CLL equation. 

A. Integrable CLL equation and soliton solution 

The well known CLL equation given as 



shares all the rich properties of an integrable system, like commuting higher conserved quantities, in- 
tegrable hierarchy, Lax pair and exact soliton solutions. Moreover due to its canonical Poisson bracket 




:■ 



a) 



b) 



iipt ~ ipxx + i2e\i/j\ 2 ip x = 0, 



(46) 
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structure this system allows an exact quantum formulation with the algebraic Bethe ansatz solution 
qDNLS. It is important to note that the CLL equation is gauge related to the DNLS equation jmp84 
and therefore it is intriguing in the present context to find the nonholonomic deformation of the CLL 
equation and examine how the gauge connection affects its relation with the deformed DNLS derived 
above. 

Since our approach for nonholonomic deformation starts with the Lax pair we derive them from those 
of the DNLS system through a gauge transformation 

Ucll(X) = h- l U dnls {\)h - h~ l h x , Vcll(X) = h^VdnMh - h~ l h t . (47) 

where the gauge field is defined as 

h = e iS(x,t)a^ with 6x = £ | 9 |2 = e |^| 2) s t = ie(^ x -^* x )+6\^\ 4 . (48) 

Using the expressions ([2]) and ([3]) we can write down explicitly the CLL Lax pair, the flatness condition 
of which would generate the CLL equation (|46p . For example, the associated space-Lax operator, which 
is crucial for the ISM application, would take the form 

Ucll(X) = ~i{\ 2 + e|Vl Vs + X(iPa+ + e^V"). (49) 

One can carry out the ISM independently for extracting the soliton solution for this equation, which 
however can be avoided by noticing a direct relation between the CLL field tp and the DNLS field q given 
by the gauge connection as 

V, = ge-W (*•*), 5(x,t) = e f \q\ 2 dy = e f \^\ 2 dy, (50) 



which constructs the CLL soliton from ([5]) and (|42h as 

. . exp(20-2icr) 

^ = 4er? (51) 

exp(4y) + exp(ze7j 

with the same expressions © for 9, a and (|43|) for the unperturbed parameters u>o and vq. The soliton 
solution under nonholonomic deformation, as we see below, can be obtained by deforming the solitonic 
parameters (|5ip in the same way as done above. 

B. Nonholonomic deformation of the CLL equation and accelerated soliton solution 

For constructing the deformed CLL equation we follow the above procedure by taking the space- 
Lax operator UcllW as (jM]), while adding a deforming part (fT4|) to the time-Lax operator VcllW- 
The deforming matrices j = 0,1,2 may be taken in the same form (]16p . though due the gauge 

transformation they are linked in this case to (|16p as — ► h~ 1 G^h. Using therefore the deformed 
Lax operators we can derive the integrable perturbed CLL equation as 

itpt ~ tpxx + 2ie\i>\ 2 ip x = 2(gip - iw), (52) 

and its complex conjugate, together with a different set of nonholonomic constraints on the perturbing 
functions: 

g x (x, t) = ie(ipw* — ip*w), w x (x,t) = 2ia(t)ip — ie\ip\ 2 w, (53) 

where w in fact is gauge related to (|19p through a transformation like ()50p . It is crucial to notice here 
that, in spite of the gauge equivalence between the unperturbed DNLS and the CLL equations, their 
respective nonholonomic deformations exhibit different structures as evident from f)53[) . As a result the 
resolution of the constraints (|53p by introducing the potential field, as done for the deformed DNLS, 
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clearly fails here, since in w x as seen from (|53p an extra term appears due to the gauge transformation. 
Therfore unlike the LF equation or the deformed DNLS discussed above, the perturbing fields can not 
be eliminated from the deformed set of the CLL equations f)52|) - f)53 [) through the local potential field 
and one has to consider (|52p as the perturbed CLL equation with nonholonomic constraints (153|) . The 
whole perturbed system however remains integrable with the associated Lax pair defined as above. The 
soliton solution for this deformed CLL equation takes again the form (|51j) . where the expressions for 
9, a should be replaced as ()38|) . Note that the unperturbed parameters are given here as (|43p . while the 
time-dependent deformed frequency uJd(t) and the soliton velocity Vd(t) are in the form (|38p . This would 
result again to an accelerating soliton, where the solution is the same as that for the deformed DNLS 
case \q\, shown in Fig. 2a. It is important to note here that in the above constraint relations (|53p . we have 
considered the same solution (|16p . which however is the simplest possible solution in the present case. 
The deforming matrix G^ in fact can have a more general structure with nontrivial off-diagonal elements 
expressed through the function e s ( x,t ' defined in (150f) . since the last constraint in (|15p is extended now 
to a more general form: G x + ze|*0| 2 [a 3 , G^] = 0. This situation would naturally make the constraint 
(|53p more complicated, yielding another novel integrable perturbation of the CLL equation. 

Exploiting the integrability of the deformed CLL equation we can construct its two-fold integrable 
hierarchy including the higher deformations, following the same recipe as above for the deformed DNLS 
equations. For this we can use again the same space-Lax operator (149)) and construct the generalized 
time-Lax operator (|12l [T3|) from those of the corresponding deformed DNLS by using the known gauge 
transformation (|47|48p . Note however that the higher order nonholonomic constraints on the deforming 
matrices GW, j = 0, 1,2, . . . are no longer same as ()10p . though may be derived from it through the 
gauge transformation (|48|) . inducing the changes like 

G U) ^ h& j) h-\ G® h(G® + [h~ l h x , G {j) })h-\ 

We however will not reproduce here such hierarchical equations explicitly. 

IX. CONCLUDING REMARKS 

The integrable nonholonomic deformation of the Kaup-Newell class including the derivative NLS 
equation, which remained as a challenging unsolved problem has been solved here completely with the 
construction of a two-fold integrable hierarchy with higher order deformations. As an important ap- 
plication we identify the recent Lenells-Fokas equation as a member of this deformed DNLS hierarchy, 
which allows us to bypass the involved procedures of LF,lenells09 and find its soliton solution as a simple 
deformation of the known DNLS soliton. The soliton solutions in such deformed systems are shown to 
exhibit unusual properties like accelerating (decelerating) solitonic motion, induced by the perturbing 
functions, which control such motions from the space-boundaries. 

The nonholonomic deformation of the Chen-Lee-Liu derivative NLS equation though can be con- 
structed similarly, gauge relation between the undeformed CLL with the Kaup-Newell DNLS equation, 
is not extended straightforwardly to their corresponding deformed equations. Nevertheless the soliton 
solution of the deformed CLL equation together with that for its deformed hierarchy can be constructed 
in a similar way. 

A promising potential application of the nonholonomic deformation for the family of DNLS equations 
arises when such deformations are interpreted as perturbations subjected to differential constraints, which 
keep the whole system integrable. Such perturbed DNLS equation can be considered therefore as a novel 
coupled system of the DNLS and a Maxwell-Bloch type equation and may have important applications 
in the fiber optics communication boosted by a doped laser system, extending the idea of NLS-MB to a 
new direction kundulO. 
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